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Abstract. We present near-infrared (NIR) surface pho- 
tometry of a sample of 14 early-type spirals with observed 
rotation curves. In this first paper, we report the results of 
two-dimensional parametric and non-parametric decom- 
positions to separate the bulge and disk components; the 
parametric bulge is modeled with a generalized exponen- 
tial law of integer index n, and the disk with a simple 
exponential. We find that the derived bulge parameters, 
for a given galaxy, vary systematically with the bulge 
shape index n. The mean early- type bulge has a best-fit 
n = 2.6, and 80% of the sample has best n of 2 or 3. Bulges 
are rarely spherical; the median bulge intrinsic ellipticity 
is 0.33. The median early-type disk has (J — K)d more 
than 0.1 mag bluer than the bulge, and a NIR disk sur- 
face brightness more than 1 mag arcsec~ 2 brighter than 
later-type disks. Our data are consistent with the well- 
established correlation of both bulge and disk surface 
brightness with physical scale length, and we note that 
the location of bulges within this projection of the fun- 
damental plane depends on their shape index n. In agree- 
ment with previous work, the ratios of bulge and disk scale 
lengths are consistent with a constant value r e /rd = 0.3; 
however, such value again depends on the bulge index n, 
implying that claims for a scale-free Hubble sequence may 
be premature. 
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1. Introduction 

Surface brightness distributions of external galaxies have 
been studied for many years but a reliable decomposi- 
tion into structural components is often difficult to ob- 
tain. The reliability of models and t echni ques has been 
questio ned by several authors ( Kent 1986 Schombert & 
Bothun 1987 ; Byun & Freeman 1995 ), and there are sev- 
eral points that need further study and refinement. First, 
brightness distributions have often been studied as one- 
dimensional (1-D) radial profiles extracted by averaging 
along elliptical annuli; such ellipses can deviate consid- 
erably from the actual isophotes, especially in highly in- 
clined systems with a luminous bulge. The result provides 
therefore a distorted profile of the surface brightness along 
the major axis. Second, the seeing must be properly taken 



into account (c.f., Schombert & Bothun 1987), especially 
when studying the central regions where brightness gradi- 
ents are highest. Third, the choice among the parametric 
forms of brightness distribution (exponentials, Hubble and 
de Vaucouleurs laws, and so on) does not depend on any 
physical argument, but only on their ability to fit the data. 
Finally, the effects of internal extinction should be taken 
into account. 

The new near-infrared (NIR) panoramic detectors have 
made imaging in the 1 to 2.5 fim regime rather straightfor- 
ward with sensitivity and accuracy comparable to those 
attainable at optical wavelengths. The NIR bandpasses 
have been advocated to be the ideal ones to study the 
characteristics of the galactic backbones, that is the stel- 
lar populations which make up the mass distribution in 
a galaxy (e.g., Rix & Rieke 1993| ). The reason for this is 
twofold: first, the extinction is lower, by a factor often be- 
tween the B and K bandpasses; and second, the emission 
of old stellar populations peaks in the NIR. 

In this context, we have undertaken a program of NIR 
imaging of bright spiral galaxies with measured rotation 
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curves. These images are used to decompose the lumi- 
nosity distribution into bulge and disk components and 
then to analyze their contribution to the observed rotation 
curves. To overcome the methodological drawbacks men- 
tioned above, we have developed a parametric technique 
to fit a two-dimensional (2-D) bulge + disk distribution 
to the entire image, which takes into account the effect 
of seeing. In addition, generalizing Kent's non-parametric 
approach, we have developed an iterative algorithm to de- 
duce from the 2-D brightness distribution the contribu- 
tions of bulge and disk, again taking into account the see- 
ing. In this first paper, we report the results of the struc- 
tural decomposition for a restricted sample of early-type 
spirals; a subsequent paper will describe the inferred mass 
distributions. Such systems appear well suited for studies 
of this kind, mainly because of the relative smoothness of 
their disks and their (reputedly) lower internal extinction. 
On the other hand, due to the low gas content, rotation 
in Sa's is not measured as far out as for later types, and 
hence they are not quite as effective for exploring dark 
matter properties. In several cases it has been noted that 
the presence of a dominant bulge might severely compli- 
cate the picture and the analysis (Kent 1988); these are 
galaxies with luminous bulges and slowly rising rotation 
curves, and some of them are also included in our sample. 

2. Observations and data reduction 

We selected 14 early-type spirals with observed rotation 
curves (Rubin et al. 1985; Fillmore et al. 1986), and with 
diameter < 8 arcmin. so that they could be imaged within 
our field of view. The sample is listed in Table |l|, ob- 
jects with anomalously slowly rising rotation curves are 
marked with an asterisk. Morphological type and coor- 
dinates are taken from NEDEl. Three galaxies reside in 
the Virg o Cluster: NGC 4419, 4450 and 4698 (Binggeli et 
al. 1985). The remaining 11 are members of loose groups 
(Huc hra & Geller|1982l Geller & Huchra|1983|; Rubin et al. 



(2.2 /im) bands, using the Arcetri NIR camera ARNICA 
(Lisi et al. |1996[), equipped with a NICMOS 3 detector. 



1985; Garcia 1993[). Distances were computed according to 



the infall model in the Third Refere nce Ca talog of Bright 
Galaxies (RC3, de Vaucouleurs et al. 199l| ), and assuming 
an H = 50kms -1 Mpc -1 . For each group or cluster its 
mean redshift was used; the adopted distance scale im- 
plies a Virgo distance modulus of 31.9 mag. The position 
angle (Column 6) is obtained at the outer isophotes, mea- 
sured East from North; the optical radius and the B-band 
magnitude (Columns 7, 8) are from the RC3; Column 9 
reports the bands in which each galaxy has been observed. 

Images were acquired at the Gornergrat Infrared Tele- 
scope (TIRGO, 1.5 m, f/20) in the J (1.2 ^m) and K 



1 The NASA/IPAC Extragalactic Database (NED) is oper- 
ated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the U.S. National Aeronautics 
and Space Administration. 



The field of view of the camera is 4x4 arcmin with a plate 
scale of 0.97 arcsec pixel -1 . Galaxies were observed by 
spending half of the total integration time (~ 15 min) 
on source, and half on the sky, alternating sky and source 
frames. The seeing FWHM was typically 2", and sky mag- 
nitudes were roughly 15.5 and 12.5 mag arcsec -2 in J and 
K, respectively. Flat fields were obtained from sky frames; 
the final images are typically flat to 0.3-0.5% in the J 
band, and to 0.1% or better in the K band, giving (la) 
limiting magnitudes of 21.5 and 20.5 mag arcsec -2 in J 
and K respectively. All image reduction was performed 
with IRAF and the STSDAS packagesB. Standard stars 
were selected from the UKIRT Faint Standard List (Casali 
& Hawarden 1992). Zero points were derived for each night 
of observations after correcting for atmospheric extinction 
using mean TIRGO values (0.11 and 0.07 m ag airmass -1 
in J and K respectively; Hunt et al. 1987). Images were 



then corrected for extinction within our Galaxy and for 
redshift. Galactic extinction was determined according to 
the Burstcin- Hcilcs values given in RC3, using the extinc- 
tion curve given in Cardelli et al. (1989); K-corrections 
for redshift were calculated from the coefficients given by 
Frogel et al. ( |1978| ). 

Color images were made by combining the J- and K- 
band magnitude images after first registering them to a 
fraction of a pixel. When the seeing differed significantly 
between the two bands, the higher-resolution image was 
degraded to the level of the lower resolution one. The K- 
band images and J — K color images are shown in the top 
panels of Fig. |l|. 

For five of the sample galaxies we can compare the 
if -band major-axis profiles with low-resolution (28" aper- 
ture) profiles obtai ned w ith a single-element photometer 
(Giovanardi & Hunt 1996 ); the data are shown as filled cir- 
cles in Fig. 0. As can be seen, the agreement is quite good 
in spite of the pointing uncertainties in the photometer 
profiles. 



3. Quantitative morphology 

3.1. The parametric decomposition 

For the 2-D parametric fits, we have modeled the sur- 
face brightness distribution with a generalized exponential 



IRAF is the Image Analysis and Reduction Facility made 
available to the astronomical community by the National Op- 
tical Astronomy Observatories, which are operated by AURA, 
Inc., under contract with the U.S. National Science Founda- 
tion. STSDAS is distributed by the Space Telescope Science 
Institute, which is operated by the Association of Universities 
for Research in Astronomy (AURA), Inc., under NASA con- 
tract NAS 5-26555. 
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Table 1. The galaxy sample 



Name 


Type 


RA. (1950) 


Dec 


(1950) 


Distance 


Pos. angle 


R25 


m B 


Band 






h 


m 


s 


O 






(Mpc) 


(deg) 


(arcsec) 


(mag) 




(1) 


(2) 




(3) 






(4) 




(5) 


(6) 


(7) 


(8) 


(9) 


N1024 


(R'jSA(r)ab 


02 


36 


30.4 


10 


37 


56 


69.2 


155 


117 


13.08 


K 


IN 2639 


(K)bA(r)a: • ayl 


08 


40 


03.0 


50 


23 


11 


70.4 


140 


55 


12.56 


J, K 


N2775 


SA(r)ab 


09 


07 


41.0 


07 


14 


35 


27.1 


155 


128 


11.03 


J,K 


N2841* 


SA(r)b: 


09 


18 


35.8 


51 


12 


31 


19.2 


147 


244 


10.09 


J,K 


N3593 


SA(s)0/a: 


11 


11 


59.2 


13 


05 


28 


17.2 


92 


157 


11.86 


J,K 


N3898* 


SA(s)ab 


11 


46 


36.1 


56 


21 


42 


34.8 


107 


131 


11.60 


J,K 


N4378 


(R)SA(s)a 


12 


22 


44.8 


05 


12 


06 


51.4 


167 


87 


12.63 


J,K 


N4419 


SB(s)a 


12 


24 


24.6 


15 


19 


24 


23.9 


133 


99 


12.08 


J,K 


N4450 


SA(s)ab 


12 


25 


58.0 


17 


21 


42 


23.9 


11 


157 


10.90 


J,K 


N4698 


SA(s)ab 


12 


45 


51.3 


08 


45 


35 


23.9 


168 


119 


11.46 


J,K 


N4845 


SA(s)ab 


12 


55 


28.1 


01 


50 


42 


27.1 


89 


150 


12.10 


J,K 


N5879* 


SA(rs)bc:? 


15 


08 


27.6 


57 


11 


25 


26.2 


2 


125 


12.22 


J,K 


N6314 


SA(s)a: sp 


17 


10 


33.1 


23 


19 


43 


139.2 


177 


43 


13.80 


J,K 


IC724 


Sa 


11 


41 


00.1 


09 


13 


10 


125.6 


60 


70 


13.40 


J,K 



bulge (Sersic |l982[ Sparks |1988| ): 
h(x,y) 

I, exp I -a v 
plus an exponential thin disk: 
Id{x,y) = /d(0)exp 




(1) 



1 



Td 



(2) 



I e and r e are effective (half-light) values, e& is the appar- 
ent bulge ellipticity, a„ is a constant relating the effective 
brightness and radius to the exponential values (see Ap- 
pendix), x and y are in arbitrary units, with origin at the 
galaxy's center, and x along the major axis. The bulge is 
assumed to be an oblate rotational ellipsoid, coaxial with 
the disk, and its apparent eccentricity e& is related to the 
intrinsic eccentricity e' b by: 



('b 



(3) 



The structural parameters have been determined by 
fitting the photometric data to the model, convolved with 
a circular Gaussian seeing disk of appropriate FWHM, us- 
ing a x 2 minimization. The fitted parameters are: the two 
surface brightnesses I e and id(0) (or /i e and /i^ when given 
in magnitudes), the two scale lengths r e and r^, the bulge 
ellipticity and the system inclination i. Because of the 
difficulty inherent in estimating a luminosity distribution 
of unknown form (bulge) which is partially embedded in 



another (disk), we did not attempt to fit the bulge index 
n explicitly (c.f. Andredakis et al. 1995| ). Instead, for each 
galaxy the bulge was modeled with four different values of 
the exponent n: 1,2, 3, and 4 (n = 4 for a de Vaucouleurs 
bulge). We then assessed the "quality" of the results for 
each set of 2-D fits, and determined the value of n that 
gave the best all-round fit, thus assigning, in effect, a shape 
index to the bulge. Three quality indicators were consid- 
ered in the assessment: the successful convergence of the 
process, the value of x 2 , and the mean absolute residual. 

When both J and K data were available, tests showed 
that bulge and disk decompositions and colors of the com- 
ponents were more stable when the two bands were fitted 
simultaneously. We also found no evidence for systematic 
scale length changes with wavelength of either component 
(c.f., Evans 1994). Hence, when two bands were available, 
we performed a single fit keeping scale lengths, inclina- 
tion, and bulge ellipticity the same for both bands, and 
letting the surface brightnesses of both components vary 
independently, but with constant color indexes (J — K)b 
and (J — K)d ■ 

The decomposition technique was extensively tested 
on a set of synthetic maps of axisymmctric bulge+disk 
distributions, covering a range of values for i, ej,, and B/D 
ratio similar to ours. We adopted an exponential thin disk, 
and generalized exponential laws for the bulge. An appro- 
priate amount of noise was added to simulate the typical 
signal-to-noise levels we achieve, and finally the models 
were convolved with a range of PSF's to simulate the ef- 
fects of seeing. The 2-D decomposition method was al- 
ways able to recover the true parameter set, within the 
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estimated errors, independently of the starting point in 
parameter space. 

We also performed a 1-D fit to the radial profile, ob- 
tained by averaging along elliptical annuli; center and po- 
sition angle of the annuli were held fixed, while the ellip- 
ticity was allowed to vary with semi-major axis. In the 
1-D case, the convolution of the model with the PSF re- 
quires i and tb to be assigned a priori: as customary in 
1-D techniques, ej, was fixed at and i was determined 
from the outer J-band isophotes. The best-fit bulge index 
n was determined in the same way as for the 2-D models. 

1- and 2-D techniques are compared in §4.1. 

The lower left panels in Fig. |l| show the results of the 

2- D decompositions in terms of cuts along the semi-major 
axis (K band). Table || gives the best-fit parameters, in- 
cluding the shape index n: the first line reports the val- 
ues from the 2-D fit, the second reports the 1-D ones, and 
the third gives the values from the non-parametric decom- 
positions (§ 3.2). Bulge and disk surface brightnesses are 
corrected to face-on values assuming optical transparency 
with C = l: 



Md = Md - 2.5 C log (cos i) 
ill = // e - 2.5 C log(l-e 6 ) 



(4) 



and are denoted by li%(K) and H%K) respectively. Ta- 
ble H lists the derived luminosities and bulge-to-disk ra- 
tios (B/D); the first line reports the 2-D parametric fit 
values, and the second gives the non-parametric values 



(see Sect. 3.2). The 1-D fixed-position- angle profiles to- 
gether with Tables || and || are available in electronic form 
on the ftp node: sisifo.arcetri.astro.it in directory 
f tp/pub/nir. 

3.2. The non-parametric decomposition 



Kent ( gglD was the first to introduce a non-parametric 
(np) method to separate bulge and disk; recently this has 
been extended to the profile analysis of 2-D images by An- 



dredakis et al. (1995). We propose here a new, completely 



2-D technique. Assuming the apparent ellipticities of the 
two components, tb and td — 1 — cos i, are known, let us 
consider the average surface brightness of an annulus of 
ellipticity, say, e^. For a sufficiently narrow annulus, A, the 
disk contribution will be Id(o,), where a is the semi- major 
axis. The total average brightness will be: 



< 1(a) > d = Id{a) + -J- / it 
b A J a 



clS 



(5) 



where Sa is the annular area. An analogous equation can 
be written for an annulus of ellipticity e^. The integral in 
Eq. (@) can be written as: 



2S A 



rir/2 



'cos 2 9 



4 sin 2 
PI 



9 d9 



(6) 



where pb and pd are the axial ratios of the components. 
We then obtain for any a, a system of integral equations: 



'Jb(a) = < I (a) > b 



2 r' 2 



Id a J cos 2 



% sin 2 9 d9 
Pd ) 



(7) 



I d {a) = < 1(a) > d 
2 r-n 



h aJ cos 2 9 + £S 



4sin 2 <9 ) d9 
Pb / 



To solve the system we first compute < 1(a) >b and 
< 1(a) >d for two sets of annuli of ellipticities e& and 
€d, respectively, each set sampling the entire light distri- 
bution. The two equations can then be solved iteratively 
for lb and I d at each a, starting with < I >b and < I >d 
as initial guesses to evaluate the integrals. 

When defining the two sets of annuli, a proper sam- 
pling of the semi-major axis is of crucial importance. A re- 
liable estimate of the integrals in Eq. [7] requires that each 
annulus of either set intersects with a sufficient number of 
annuli of the other set. If we constrain this number to be 
the same for every ellipse, the sampling is linear in log a, 
i.e. the annuli become wider with increasing radius. 

Results of the np decomposition technique are given 
in the third line of Table |2|, and in the second line of 
Table |[ The non-parametric parameters are defined as 
follows. The bulge r e is the semimajor axis of the ellipse 
enclosing half light; fi e is the bulge surface brightness at 
r e ; rd is the semimajor axis of the ellipse enclosing 0.264 
of the disk light, as in the exponential case. Since the 
central disk is not tightly constrained by the np fits, p,d is 
the central brightness of an exponential disk having scale 
length rd and the same total luminosity as the np disk. 
The np colors of the components are derived by the total 
component luminosities in the two bands. 

3.2.1. Testing of non-parametric decomposition 

The np decomposition algorithm was tested on the same 
set of synthetic maps used for the parametric tests. Since 
the method applies only to components with different ap- 
parent ellipticity, we only tested models with eb at least 
10% smaller than ed- Under these conditions, to within 
the noise level, the method always proved to be able to 
recover the true bulge and disk distributions given the 
correct values of the ellipticities. 

When errors on the ellipticities are introduced, the 
shape of the resulting components is affected in a sys- 
tematic way: the value of ej, affects mainly the inner part 
of the disk, while td determines the shape of the outer 
bulge. In particular, if tb is underestimated the inner disk 
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Fig. 1. if-band images, J — K color images, parametric, and non-parametric decompositions. The upper left panel shows the 
iT-band map; contours range from 21 to 14 mag arcsec - in unit steps. Images are oriented with North up and East left, and 
offsets are given in arcsec. The upper right panel shows the J — K color map, with the color table ranging from 0.5 (black) 
to 1.3 (white). The lower left panels show the results of the parametric decomposition with the K-hanA elliptically-averaged 
radial profiles shown as filled triangles, and the major-axis cut shown as a solid line. The parametric bulge and disk are shown 
as dot-dashed and dashed lines, respectively, and their sum as a dotted line. 28" aperture profiles (Giovanardi & Hunt 1996) 
are shown as large filled circles, if available. T he rem aining lower left panels show the elliptically-averaged profiles in J — K, 
and r — K with the r-band taken from Kent ( 1988| ), if available. Also shown in the J — K plots are the fitted color profile 
(dotted line), and the fitted colors of bulge (dot-dashed line) and disk (dashed line). The lower right panel shows the results of 
the non-parametric decompositions, with the same coding as for the parametric decompositions. 



Table 2. Decomposition results: structural parameters 



is too steep, whereas if the bulge is too elliptical the inner 
disk develops a hole which gets deeper as ef, increases. If 
€d is too small the outer bulge is too steep and viceversa. 
Although qualitatively these are expected trends, they are 
difficult to quantify in practice, since they depend on the 
shapes of the distributions, on the cllipticitics, and on the 
seeing. Especially when the two components have similar 
ellipticities, even small errors on t h and n d (< 0.1) per- 
ceptibly affect parameters such as the luminosity of the 
components or the B/D ratio, as well as the shape of the 
profiles. 

3.2.2. Seeing corrections 

Simulations show that, if no correction is applied, see- 
ing significantly affects np decompositions. The effects are 
stronger in the inner parts, the main feature being a cen- 
tral hole in the disk, whose extension is proportional to 
the seeing width. In fact, the seeing makes inner isophotes 
rounder, causing the component of higher ellipticity to be 
depressed. The disk profile, however, can be corrected by 
extrapolating the outer points to the center. 

To recover the intrinsic bulge profile, we introduced a 
second iterative algorithm which proved to be quite effec- 
tive on synthetic maps. The effect of seeing is accounted 
for by defining a correction coefficient fcf,(a), for each point 
of the bulge profile, as the ratio of the convolved bulge 
distribution and an estimate of the true one. Each decom- 
position is performed using fcb(a) to degrade the assumed 
true distribution. After decomposition, the coefficients are 
recomputed by convolving the extracted bulge profile with 
the PSF, and a new decomposition is performed until both 



the coefficients and the profiles converge. The initial esti- 
mate of k{,(a) is provided by the convolution of the bulge 
profile extracted without any correction. 

3.2.3. Estimate of the ellipticities 

To determine the apparent ellipticity of the components, 
we extracted from each image a radial profile of the ellip- 
ticity by fitting isophotes of variable ellipticity and fixed 
position angle. We then selected the lowest value beyond 
one seeing disk from the center as ef,, and the average value 
in the outer parts as €d (see also Andredakis et al. 1995). 
The values for the J and K bands were then averaged 
to obtain a single value of ef, and td for each galaxy. The 
accuracy of these estimates is difficult to assess: bright 
bulges can affect the isophotes even at large radii, lead- 
ing to an underestimate of e^, while bars and spiral arms 
may often lead to overestimates. In turn, ef, can be easily 
underestimated because of the seeing. From simulations, 
it seems that a reasonable estimate for the uncertainty in 
e<2 is about 0.1; it is certainly worse in ef,, 0.2 or more in 
the worst cases. 

We note that comparing the observed image to the ex- 
tracted bulge+disk distribution, for example in terms of 
residuals, is of little use to test the ellipticities, since it is 
often possible to find excellent agreement for a wide range 
of ef, and e^. A more stringent constraint is provided by 
the color profiles of bulge and disk, which turn out to be 
quite sensitive to the values of ef, and e^. This has ulti- 
mately allowed us to verify the ellipticities a posteriori (in 
particular ef,) by letting them vary within the estimated 
errors and checking the plausibility of the resulting color 
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profiles. The final accuracy is estimated to be better than 
0.05. 



Table 3. Decomposition results: _K"-band integrated luminosi- 
ties and bulge-to-disk ratios 



are not significantly different from np ones: the mean ra- 
tio of the np eh's and those fitted by the 2-D parametric 
method is 1.1 ± 0.5. System inclinations agree very well, 
with a mean ratio of 0.97 ± 0.09. Also, the bulge ellip- 
ticities determined by Kent (1988) are about the same as 
those used in our np decomposition, but with large scat- 
ter; the inclinations are in good agreement, with a mean 
ratio of 1.01 ± 0.09. 



4-2. Dependence on n of parametric components 



4. Results 

4--1- Accuracy and comparison of the different decomposi- 
tion methods 

As Fig. [j] clearly shows, the elliptically-averaged profiles 
reproduce quite well the major-axis cut in almost all cases. 
In fact, our sample contains only moderately inclined sys- 
tems (i < 75°), whose isophotes are reasonably elliptical. 
However, a comparison between the 1-D and 2-D results 
from Table || shows that the 1-D decomposition is always 
less accurate, and sometimes inconsistent with the 2-D 
one. The inferior quality of the 1-D fits is attested to both 
by the larger errors on the parameters and by the larger 
values of x 2 - A comparison between 1-D and 2-D results 
for the galaxies with significant non-axisymmetric struc- 
ture suggests that 2-D fits manage to reproduce the whole 
distribution, whereas 1-D decompositions tend to match 
closely only the inner part of the profile, i.e. the less noisy 
one (see also Byun & Freeman |1995 ). Because of their 
superiority, when discussing parametric models, we will 
consider only the results of 2-D decompositions. 

No systematic differences are found between the vari- 
ous fitted parameters: the 1-D bulge and disk scale lengths 
are, to within the scatter, the same as those of the 2-D fits, 
and the fitted inclination of the 2-D model corresponds, 
on average, to the fixed value determined from the outer 
isophotes (and used in the 1-D fits). However, disks are 
definitely better constrained than bulges: the rms differ- 
ence of the 1- and 2-D fid's (0.6 K-m&g) is half that for 
fi e (1.1 -KT-mag), and the dispersion in the ratio of 1- and 
2-D r e 's is 80%, while that for the r^'s is four times lower 
(21%), comparable to the discrepancies among disk scale 
lengths from different authors found by Knapen & van der 



Kruit p9l| ). 

We have compared the results of our non-parametric 
(np) decompositions with the 2-D parametric ones, and 
with those of Kent (1988). There is no systematic differ- 
ence between our 2-D and np bulge and disk luminosities, 
as the mean difference is —0.4 ± 0.6 mag for Mb, and 
—0.1 ± 0.4 for the disk. Again, disk parameters are more 
consistent than those of the bulge. Parametric ellipticities 



In order to compare our results with those of other au- 
thors, we have investigated trends with n of the bulge 
and disk fitted parameters. Figure ^| shows the systematic 
variation in these, for a given galaxy, as a function of n. 
It can be seen from the figure that the derived bulge pa- 
rameters strongly depend on the form of the fitting func- 
tion. The same bulge, when fitted with small n, appears 
to be "denser" (that is to say brighter fi e ), more compact 
(smaller r e ), and less luminous than when fitted with large 
n. Quantitatively, the changes are dramatic as n goes from 
1 to 4 with a mean change in fi e of 3 mag arcsec -2 , and in 
scale length of roughly a factor of 3. Moreover, the disper- 
sion in the fitted values also increases with n; the spread 
of log r e at n = 4 is 1.5 times larger than that at n = 1, 
and the spread of fi e is more than twice as large. 

The derived disk parameters also change with the n of 
the bulge, fid tends to be fainter for bulge n larger, and, 
as for the bulge, the dispersion in fid increases with n. 
The only parameter that is stable with n is the disk scale 
length rd, although its dispersion does increase slightly 
with n. 

We conclude that, at least statistically, bulge struc- 
tural parameters are strongly influenced by the form of the 
function used to derive them. Not only are the parameters 
themselves altered by constraining the form of the bulge, 
but also the dispersion in the parameters is changed. In- 
dependently of the best- fit n, requiring a de Vaucouleurs 
law to fit the bulge yields more tenuous, extended, and lu- 
minous spheroids, together with wider distributions of the 
parameters (larger dispersion), than does using a simple 
exponential. 



4-3. Bulges 



In terms of the "quality indicators" mentioned in Sect. 3.1 
n = 3 gave superior results for the majority of galaxies, 
while n = 2 was the second-best choice. In only two cases 
did n = 4 give the highest quality fit, and in one case 
(NGC 3593) n = 1. These values are in general agre ement 
with the trend noted by Andredakis et al. ( 1995 ) who 
fitted the n of non-parametric bulge profiles and found 
that early- type spirals tend to have bulges with n ~ 2 — 3. 
We also note that, in our experience, it is usually difficult 



Moriondo et al.: NIR photometry of early-type spirals. I 



7 



Fig. 2. Bulge and disk fitted parameter variation with n. Bulge parameters are shown in the left panels: from top to bottom, 
absolute K- magnitude Mb(K), effective surface brightness nl(K), and effective scale length r e . Analogous disk parameters are 
shown in the right panels. The curves show, for a given galaxy, fit parameters obtained with n going from f to 4. H%(K) and 
H C d(K) are corrected to face-on and in units of mag arcsec -2 . Scale lengths are in arcsec. 



to determine the best n with a precision much better than 
1. 

The distribution of the best-n bulge parameters is il- 
lustrated in the left panels of Fig. 0. The median bulge 
parameters are H%{K) = 16.8 mag arcsec -2 , and r e = 1.6 
kpc, as reported in Table ||. Our bulges are more tenuous, 
larger, and more luminous than those of similar type in 
de Jong ( 1996a| ), as expected given the fixed n — 1 used 
by him. The median apparent e& is 0.24 from parametric 
decompositions and 0.25 in np decompositions. This trans- 
lates into a median intrinsic ellipticity of 0.36 or 0.33 from 
np values. Notably, 0.33 is also the commonest intrinsic el- 
lipticity found in elliptical galaxies, if they are assumed to 
be rotational ellipsoids, either oblate or prolate (Mihalas 
& Binney 1981). In any case, bulges are rarely spherical 
and the results of studies assuming so should be treated 
with caution. 



Table 4. Statistics of parametric components 







Mean (a) 


Median 


Ran 


ge 




n 


2.6 (0.8) 


3 


1 - 


4 






1.9 (1.6) 


1.6 


0.4 - 


6.4 


Bulge 




16.7 (1.1) 


16.8 


15.1 - 


18.9 




tb 


0.24 (0.13) 


0.23 


0.04 - 


0.50 




{J-K) b 


1.13 (0.27) 


1.06 


0.72 - 


1.55 




M b (K) 


-23.9 (1.0) 


-23.8 


-25.3 - 


-22.2 






4.9 (2.8) 


4.6 


1.40 - 


11.4 


Disk 




17.3 (0.78) 


17.1 


16.1 - 


18.9 




(J~K) d 


0.96 (0.08) 


0.94 


0.85 - 


1.11 




M d {K) 


-24.4 (0.8) 


-24.3 


-25.4 - 


-22.5 



Correlations between (average) surface brightnes s and 
scale length ha ve be en found for spiral bulges (Kent 
Kodaira et al. 1986 ) and for ellipticals (Korm cndy 
Hoessel 



1985 



1977 



Schneider 1985 ; Djorgovski & Davis 1987 ), and 
these two observables constitute an almost face-on view 
of the "fund amental plane" (FP) (e.g., Kormendy & Djor- 
govski 1989] , and references therein). The left panels of Fig. 
H show scatter plots of bulge ^ c e {K) and < ^ C {K) > e vs. 



r e , < fJ. c (K) > e is the average surface brightness within 
the half-light isophote commonly used in FP studies. The 
upper panel shows results for all values of n, and the lower 
one only best-n values. Our results are consistent with the 
s lope within the FP found for bulges by Andredakis et al. 
( 1995 ), shown as a dashed line in the figure. Although 
we have transformed their regression line to our distance 
scale and to the i\-band according to their precepts, we 
note a slight offset: our best-n bulges are generally dimmer 
(~ 0.5 mag), for a given r e , than theirs. 

Figure || gives further insights as to what happens 
when the form of the bulge is constrained to one value 
of n. As discussed in the previous section, when for ex- 
ample n is 4, the resulting bulge parameters lie in the 
tenuous, extended portion of the parameter space; when 
n is 1, resulting bulges are smaller and denser (see also 
Fig. |^). From Fig. [|, it appears that different n values, 
that is to say different bulge shapes, occupy different re- 
gions of the FP, such a behavior is evident in both the 
upper- and lower-left graphs. An analysis of the correla- 
tion coefficients shows that while the best-n set of points 
is significantly correlated, as is the global set of points for 
all values of n, each individual set with fixed n is not. 
Moreover, the slopes of each fixed n group increase with 
n, although even the slope of n — 4 is not as large as the 
global one. The appearance of the top-left graph is deter- 
mined mainly by a "geometrical" effect, that is a constant 
luminosity relation although, see Fig. ||, higher n's pro- 
duce slightly more luminous bulges. Independently of the 
details, it is clear that the large scatter related to the un- 
certainty on the decomposition has a high incidence on 
the position of a bulge within the FP, as shown by the 
error bar in the lower-left panel. 

It has been suggested that residuals relative to the FP 
are correlated with shape parameters (Hjorth & Madsen 
1995; Prugniel & Simien 1997). Although this projection 
of the FP is not appropriate for such considerations, the 
lower-left graph (bulge best n) in Fig. |i| suggests that even 
the distribution within the FP is at least partially gener- 
ated by form variations. The distribution of the bulge np 
parameters does not reveal any dramatically different be- 
havior; if anything, the distribution is tighter and situated 
in the low-n region of the plot. 
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Fig. 3. Distributions of best-n bulge and disk parameters. Surface brightnesses are corrected to face-on and in units of 
mag arcsec" . Scale lengths are in kpc. 



Fig. 4. Scatter plots of bulge and disk ^ C (K) vs. logr. Scale lengths are in kpc. The left panels show bulge parameters, and 
the right those of disks. The upper panels show points for all values of n, with n = 1, 2, 3, 4 shown by x's, open circles, filled 
triangles, and filled squares, respectively. The lower panels show only the best-n points, with n denoted as in the upper panel, 
together with the non-parametric results shown as crosses. In the lower panels, < /J. C (K) > e is the average surface brightness 
within the effective isophote ; note that also for disks the abscissa is r e . In the lower panels, the dashed line shows the regression 
found by Andredakis et al. (1995) converted to our distance scale and to the A"— band according to their precepts. The inclined 
error bars represent the average shift related to an uncertainty of ±1 in the bulge index n. 



4.3.1. Bulge colors 

The median (J — K)b of 1.06 (1.04 for the np decom- 
position) agrees with the colors measured by Giovanardi 
& Hunt ( |1996| ), and is redder by about 0.1-0 .2 mag than 
those measured in later types (Frogel 1985; Giovanardi 
& Hunt 1988). The scatter is large, 0.3 mag, with some 
bulges having J - K as high as 1.5 (NGC 4845). We find 
that (J — K)b correlates (98% significance) with ^(K), in 
the sense that redder colors are associated with "denser" 
bulges. In contrast, (J — K)\, is independent of Mb(K), 
and of total galaxy luminosity. 

The four objects (NGC 3593, 4419, 4845, and IC 724) 
with ( J—K)b > 1.3 also show red extended circumnuclear 
structure in the color images of Fig. |l|, inflections or bumps 
in their surface brightness profiles, and red gradients in 
the inner color profiles. Such features have been observed 
in starburst galaxies (Hunt et al. 1997), and we would 



argue that the red J — K bulge color is revealing star for- 
mation in progress. The most clear-cut case is NGC 3593 
which, besides a h igh m id-infrared 12 /zm surface bright- 



ness ( Soifer et al. 1989 ) and high molecular gas content 
(Sage 1992 ), hosts two counterrota ting s tellar disks and 



a disk of ionized gas (Bertola et al. 1996 ). Moreove r, Ha 
images reveal an HII- region ring (Pogge & Eskridge 1993 ) 
whose structure closely resembles that seen in our J — K 
image. NGC 4419, the only barre d galaxy in our sample, is 
a LINER (Huchra fc Bu rg 1992 ) with m id-infrared prop- 
erties (Soifer et al. |1989 Devereux 1987 ) and CO content 
(Young et al. 1995 ) typical of starburst s. NG C 4845 was 
defined as a starburst by David et al. ( 1992 ) on the ba- 
sis of its FIR-to-blue luminosity ratio and X-ray excess. 
IC 724, one of the most distant in our sample, harbors 
more than 10 9 M Q of HI (Eder et al. 1991), but we have 



found no evidence in the literature for star formation ac- 
tivity. The J — K image may be just revealing a normal 
bulge, partially obscured by a dusty disk. 



4-- 4- Disks 



The distribution of the disk parameters is shown in the 
r ight p anels of Fig. 0. Similar to the results of de Jong 
( 1996a ) for early spiral types, the median disk has a fJ.^(K) 
of 17.1 mag arcsec^ 2 , r d = 4.6 kpc; with M d (K) = -24.3 
mag it is slightly more luminous than the median bulge. 
The median ratio of and isophotal (optical) radius R25 
is 0.24 (shown in Fig. | as a dotted line in the upper 
left panel), comparable to what is found in late- type spi- 
rals (Giovanardi & Hunt |1988| , Giovanelli et al. |1995| ). Al- 
though similar in size, these early-type disks are more than 
1 K-va&g arcsec -2 brighter than those in late-type spirals 
(Giovanardi & Hunt |1988| ). 

The right panels of Fig. |] show scatter plots of disk 
Hd(K) vs. rd (upper panel) and of disk < fJ. c (K) > e vs. r e 
(lower panel). As for the bulge, correlations of disk u£(K) 
with rd have been noted for some time (e.g., Kent 1985 ). 
It is interesting to note that, when plotted in terms of 
the photometric observables commonly used in FP stud- 
ies (lower-right panel) , disks dwell in a region of this FP 
projection which is contiguous and similar in shape and 
extent to that of bulges. The disks appear to extend the 
bulge relation to larger radii and fainter surface bright- 
nesses. Also evident, in the lower-right panel, is the rough 
consis tency with the slope for bulges found by Andredakis 
et al. ( 1995 ), although with a large offset. It is clear from 
Fig. |^ that, unlike the bulge, the relation between disk 
parameters does not vary substantially with bulge n. 
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4.4.1. Disk colors 

The median (J — K)d of 0.94 (0.91 for the np decompo- 
sition) is similar to the central colors of late-type spirals 
(Frogel p85| ; Giovanardi & Hunt |1988[ ). The scatter about 
the mean is 0.08 mag, smaller than for (J — K)^. The me- 
dian disk is 0.12 mag bluer than the bulge, an effect not 
noted by Terndrup et al. (1994) whose sample was domi- 
nated by later types. As for bulges, redder disks tend to be 
"denser" (98% significance ), and (J — K)d is independent 
of M c ] and inclination. 



4.4.2. Opacity of the disks 

Although the diagnostics of dust content in galaxies have 
been extensively revised in recent years (e.g. Byun et al. 
1994 - BFK; Bianchi et al. 1996| ), such studies have made 
clear that disk opacities are not easy to determine. In 



the following we gather the indications about disk opacity 
obtained here; none of the tests we adopt is particularly 
stringent, due mainly to the small number statistics, but 
all converge on conservative estimates for the opacity of 
early type disks: ry(0) ranges from 2 to 4, where Ty(0) is 
the central optical depth in the V band (face-on). This is 
essent ially the same result reach ed by Peletier & Willner 
(|1992|) and Giovanardi & Hunt ( |l996|) . 

(a) Correlation between apparent disk brightness and 
inclination. We find a slight trend in both bands, with 
slopes Cj = 0.66 ± 0.23 and C K = 0.73 ± 0.23, both 
compatible with a fully transparent disk (with C = 1, 
see Eq. 4). Taken at face value, a C — 0.7 corresponds 
to a Ty(0) of ~ L.1 if measured in the J band, and to 
Ty(0) ~ 1.8 if in These moderate values for the central 
opacity imply that the spread observed in the NIR \i c d 
is intrinsic and not due to extinction; de Jong (1996b) 
reached a similar solution on the basis of a larger sample. 

(b) Variation of disk scale length with wavelength and 
inclination (Evans 1994 ; Peletier et al. 1994 ). Five of the 
sample galaxies have been paramctrically decomposed in 
the o ptical : either r (Kent |1985| NGC 2639), V (K odaira 
et al. |1986|; NGC 3898, 4698), or B (Boroson gggl NGC 



2775, 2841, 3898). We find a trend in the ratios of our to 
their r^: for the only measurement in r the ratio is exactly 
1, it decreases to 0.85 in V, and to 0.70 in B. In addition, 
these ratios depend on the inclination, thus providing a 
test which is largely free from the influence of intrinsic 
color gradients. The correlation, in the sense of smaller 
ratios for higher i, implies ry(0 ) < 3 . These results are 
consistent with Peletier et al. ( |1994|) , who find that B 
and K scale length ratios vary from 1.2 to 2.0, and with 
inclination. 



(c) Colors. As noted in the previous section, (J — K)d 
does not depend on i. We estimate the maximum (3 a) 
slope of (J — K)d vs. sec i which is still compatible with 
our data to be 0.075. For a Triplex model (Disney et al. 
19890 with C = 0.5 (Peletier & Willne r |l992T ) , this implies 
a Ty(0) < 3. We noted in Sect. 4.4.1 that red disks were 



associated with bright fi d , which again points to moderate 
opacities. Indeed, since \i c d is corrected for inclination as- 
suming transparency, a high opacity would translate into 
faint brightnesses for reddened disks. 

4-- 5. Relationship between bulge and disk 

The best-n median B/D ratio is 0.8 with values ranging 
from 0.2 to 2; even in this early-type sample, more than 
two thirds of the galaxies have disks more luminous than 
bulges. With the exception of NGC 1024, B/D ratios ob- 
tained from the np decomposition are always less than 
1 (as can be seen from the lower left panel in Fig. [|). 
The two methods yield B/D values which differ by al- 
most a factor of ~ 2, but with large scatter. We have 
verified that this is mainly imputable to the choice of ef, 
and i, the values adopted in the np case being lower. Our 
parametric B/D ratios are comparable, although some- 
what larger, to those (parametric) found by Kent ( 1985 ) 
in the r band. Also, our B/D's (both parametric and np) 
in the K band are 10 ~ 15% larger than in J. That the 
if -band B/ D is la rger than in the optical was also noted 
by de Jong ( 1996a ), but with values smaller than ours due 
to his choice of n = 1 bulges. It is evident that the B/D 
ratio is rather model dependent, and different decompo- 
sition methods provide estimates differing by factors of 2 
or more, as illustrated in Fig. [s] where logr e , log r e /rd, 
and log BLD are shown as a function of optical (isopho- 
tal) radiusQ. Inspection of this figure also shows that the 
derived bulge and disk parameters, including best n, are 
not appreciably affected by biases associated with galaxy 
apparent size. 

A linear correlation between r e and ov er all spiral 
types has been recen tly found by de Jong ( 1996a ) and 
Courteau et al. ( 1996] ). They interpret the correlation as 
an indication that the Hubble sequence is scale-free since 
the relative size of bulge and disk does not depend on 
morphological type. It can be seen from the middle pan- 
els that our data are also consistent with constant r e jvd\ 
the sample median best-n r e /rd of 0.3 is shown as a dot- 
ted line. This value is a factor of 2 larger than that found 
by Courteau et al. in the r band, (r e /rd) = 0.13j-.and by 
de Jong for the K' data alone, (r e /rd) = 0.14 a. How- 
ever, r e /rd appears to be strongly influenced by the bulge 



3 Fro m Fig . 8 in BFK and assuming A(X) /A(V) as in Cardelli 
et al. ( 1989T) ; these estimates have been corrected approxi- 



4 These were taken from Table 

5 Their values have been converted to bulge effective scale 



mately for the lower dust albedo in the NIR. 



lengths using the conversion factor for n 
Appendix. 



1 given in the 
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Fig. 5. Disk scale lengths, ratio of bulge and disk scale lengths, and B/D vs. optical isophotal radius i?25- Best-n values are 
shown in the left panels, and values for all n are shown on the right; symbols are as in Fig. The upper panels show logr^ vs 
log7?25; the dotted line in the upper left panel illustrates a linear relationship between ra and -R25 with a mean ratio rd/R.25 of 
0.24. The middle panels show the ratio \ogr e /rd vs logi?2s; the dotted line on the left gives the best-n median constant value 
Tc/fd = 0.3, while the pair of lines on the right show the n = 1 median r e /rd = 0.2, and the n = 4 median r e /rd — 0.7. The 
lower panels show log(£?/D) vs logi?25- 



parameterization: the dashed lines shown in the middle 
right panel of Fig. || illustrate the values obtained from our 
n = 1 fits (sample median r e /r c i — 0.2), and for the n — 4 
fits (median r e /rd = 0.7). According to whether bulges are 
fit with a simple exponential or with the de Vaucouleurs 
law, r e /rd changes by more than a factor of 3. Hence, if the 
best-fit n changes with m orphological type as suggested by 
Andredakis et al. ( 1995 ), the claims made by Courteau et 
al. for a scale-free Hubble sequence may be premature. 

4-6. Color gradients 

Since we adopt the same scale length in J and K, our 
parametric decompositions yield bulges and disks with 
uniform color. On average, the resulting bulges are red- 
der than the disks by more than 0.1 mag, and we should 
detect a significant color gradient at the transition be- 
tween bulge and disk. Such gradients are clearly evident 
in NGC 3593, 3898, 4419, 4845, 6314 and IC 724, all ob- 
jects whose bulge and disk colors differ greatly. When such 
gradients are present, they also appear, enhanced, in the 
r — K profiles. 

Regarding the gradients within the single components, 
we give no estimate of bulge color gradients! For the disk, 
following Terndrup et al. ( [1994 ), we estimated outer (> 3 
kpc) color gradients, computed by fitting J — K and r — K 
versus logr (in kpc); thejtjwill be denoted with 8{J — K) 
and S(r — K) respectively!! We detect J — K gradients at 
3<7 in only one case, NGC 6314, and in six if we consider 
a 2 a limit. All galaxies with J — K gradients for which 
we have r data, namely NGC 3593, 3898, 4378, 4419 and 
6314, also show significant S(r — K). Significant 8(r — K) is 
also found in NGC 26 39, 4845 and IC 724. In agreement 
with de Jong ( 1996b ), we only find negative gradients, 



6 This could be accomplished, in principle, using np profiles 
or the inner part of the obseved color profiles but they are 
difficult to detect. Due to the seeing, the very inner color pro- 
file is hardly reliable, while at outer radii the bulge brightness 
rapidly falls deep below the disk. 

7 We used global profiles rather than np disk profiles because 
the latter are not available for r — K, and because the two types 
of profile are very similar in J — K. 



ranging from —0.24 to —0.47 mag per decade in J — K, 
and from —0.26 to —1.03 in r — K. There is no correlation 
between 5(J — K) and inclination, but inclined galaxies 
have steeper r — K profiles (see Fig. ||). 

Such color gradients provide a last assessment of the 
disk opacity, already discussed in § 4.4.2 . While NIR colors 



are stable across the disk, we find a prevalence of negative 
trends in r — K, especially in inclined galaxies (see Fig. [3]). 
A weighted fit of S(r — K) vs. —2.5 log (cos i) yields a 
slope of —0.30 ± 0.08. In B — I, a common feature of the 
models (e.g. BFK) is that gradients tend to steepen with 
increasing i only in rather transparent disks; for ry (0) > 3 
the trend is reversed due to a saturation effect. It has 
also been shown by Bianchi et al. ( |1996| ) that disk color 
gradients are not greatly influenced by the dust scattering 
properties. In r—K the reddening will be larger by a factor 
of 2 for the same 7y(0), so the observed slope is roughly 
indicative of a ry(O) ~ 1.5 (BFK; Bianchi [1995^. 



Fig. 6. Colors and color gradients vs. inclination. The upper 
panels show the mean of J — K and r — K colors beyond 3 kpc 
from the center. Color gradients evaluated in the same region 
are shown in the lower panels. 



5. Summary 

We have decomposed J— and AT— band images of 14 early- 
type spirals into bulge and disk components. 2D non- 
parametric solutions and results from fitting a parametric 
model of generalized exponential 1/n bulges and simple 
exponential disks are compared, and general characteris- 
tics of early- type spiral bulges and disks are examined. We 
find that: 

1. Even using objective and refined techniques, the de- 
composition in structural components is far from being 
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a robust and unique process. For the parametric meth- 
ods, significantly different decompositions are obtained 
for different bulge distribution laws. Non-parametric 
techniques, on the other hand, appear to be affected by 
the choice of the ellipticities of the components which 
are difficult to evaluate objectively. Of the two compo- 
nents, the bulge is the most subject to errors, since the 
inner part is masked by seeing, and the outer regions 
are buried beneath the disk. 

2. Bulge structural parameters are strongly influenced by 
the form of the function used to derive them. The same 
bulge, when fitted with small n, appears to be "denser" 
(brighter fi e ), more compact (smaller r e ), and less lu- 
minous than when fitted with large n. The dispersion 
of the fitted parameters also increases with n. 

3. The median early- type bulge has a shape index n be- 
tween 2 and 3, nl(K)=16.8 mag arcsec -2 , and r e =1.6 
kpc. It is also red, with (J — K)b — 1.06, and redder 
bulges tend to be "denser", that is with brighter /x e . 

4. As noted by Kent ( 198S| ), bulges are rarely spherical. 
The median intrinsic ellipticity is 0.34, equivalent to 
a disk with 50° declination. This restricts the applica- 
bility of non-parametric techniques to rather inclined 
systems and suggests treating with caution the studies 
which assume spherical bulges. 

5. The median early-type disk with pt d {K) = 17.1 is more 
than 1 mag arcsec" 2 brighter than later-type disks, 
and bluer than the bulge in (J — K) by more than 
0.1 mag. Disk scale lengths agree fairly well with those 
found by other at different wavelengths, and we con- 
firm a tendency for NIR disk scale lengths to be smaller 
than those at optical wavelengths (e.g., Peletier et al. 

Td/R2b is approximately constant, 0.24, simi- 



G 



1994|) 



lar to the value of 0. 25 fou nd for late- type spiral disks 
(Giovanardi & Hunt 1988 Giovanelli et al. 



19951) 



Both bulge and disk surface brightnesses correlate with 
respective scale lengths, consistently with the projec- 
tion of the fundamental plane fo r ellip ticals and spi- 
ral bulges (e.g., Andredakis et al. 1995[ ). We note that 
uncertainties in the decomposition, especially in the 
shape index n, strongly influence the position of a 
bulge within the FP. Disks appear to reside in a re- 
gion of this FP projection which is roughly contiguous 
to that of bulges, extending the correlation to larger 
radii and fainter surface brightnesses. 
7. We confirm the tendency for the ratio of bulge and 
disk s cale le ngths r e /ra to be co nstant , noted by de 
Jong (1996a) and Courteau et al. ( |1996| ). However, we 
find a mean (best-n) value r e /rd = 0.3, significantly 
larger than the value found by de Jong and Courteau 
et al.; our n = 1 value of r e /rd of 0.2 agrees roughly 
with their value of 0.13-0.14, while our n = 4 value 
is 0.7, more than a factor of 3 larger. We attribute 
such differences to different bulge parameterizations 



and caution that if best n varies with morphological 
type, as suggested by Andredakis et al. ( 1995 ), r e /ra 
may not be constant with morphological type, and thus 
the Hubble sequence may not be scale free as proposed 
by Courteau et al. (1996). 
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A. Generalized exponential brightness distribu- 
tions and conversion factors 

Assuming circular symmetry, we define a generalized ex- 
ponential bulge as a surface brightness distribution: 



J(r) = Jq exp 



l/n 



(Al) 



where r is the projected radius, Iq the central brightness, 
ro the exponential folding length, and n an integer index 
(n > 1). With the change of variable x n = r/ro, it is easily 
seen that the total luminosity L is: 



L = 



I(r) 2nr dr = n(2n)\ Vo 



(A2) 



It is customary to introduce instead of ro an effective 
radius r e , which encircles half of the luminosity, and in- 
stead of Io an effective brightness I e = I {re)- We derive 
in the following the relations between these quantities in 
the case of a generalized exponential. 

The relation defining r e is: 



I(r) 2nr dr = — 



(A3) 



with the above mentioned change of variable, this integral 
yields: 



«7 



2n— 1 

^ (2n 



1 



i=0 

which we write as: 



= — , with 



- , (A4) 



m a 
i=0 



with to = 2n .— 1. 



(A5) 



The left hand side of this last equation is a cumulative 
Poisson distribution with parameter a n . This can be writ- 
ten in terms of x 2 probability functions (Abramowitz & 
Stegun |l97l|) : 



i=0 



(A6) 
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with x 2 = 2a„ and v — 2(m + 1) = in. 

a n is therefore easily evaluated from \ 2 tables. In any case, 

for large v values, (x 2 /^) 1 ^ 3 is approximately normally 

distributed: 



Cinzano P., Corsini E.M., et al., 1996, ApJ 458, 



Q(x 2 



with X2 



Q{X2) , 



(A7) 



(A8) 



Since in our case Q{x2) = \, it follows x-i ~ and \ 2 
v(\ — j^) . Then, with the approximations adopted: 



ro ( 2n - - 



(A9) 



It turns out that the approximation is already quite good 
for n = 1, with a relative error of 0.7%. Table || reports 
the solution of Eq. AE and its approximation with Eq. A9 
for various values of n. 

For the effective surface brightness we have: 



J(r e ) = Io exp 



I e a " ~ / exp ( — 2n + 



(A10) 



These results are easily generalized to elliptical distribu- 
tions and they are found to hold in the same form if r e is 
the major semiaxis of the ellipse encircling half luminosity. 



Table 5. r e approximation 



11 


V 


On 


2n-| 


1 


4 


1.6784 


1.6667 


2 


8 


3.6721 


3.6667 


3 


12 


5.6702 


5.6667 


4 


16 


7.6693 


7.6667 


6 


24 


11.6684 


11.6667 


10 


40 


19.6677 


19.6667 
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Table 2. Decomposition results: structural parameters 



Name 




(J-K) e 


r 


e 


e-b 




(J - K) d 




i 




n 






(2) 


(3) 


(4) 


(5) 


If. 

y° 




(7) 


(8) 


(9) 


(10) 


(11) 


NGC 


1024 


17.33 


(0.04) 




3.16 


(0.13) 


0.36 (0.01) 


18.85 


(0.14) 




8.21 (0.49) 


66.7 ( 1.2) 


1.2 


3 






17.28 


(0.11) 




4.79 


(0.40) 










7.02 (0.01) 


68.9 


26.5 


















0.35 










69.3 


1.0 




NGC 


2639 


18.84 


(0.07) 


0.84 


6.42 


(0.18) 


0.32 (0.01) 


16.11 


(0.06) 


0.98 


3.38 (0.03) 


65.2 ( 0.3) 


1.2 


4 






16.59 


(0.46) 


0.82 


1.55 


(0.60) 




15.71 


(0.13) 


0.94 


3.49 (0.16) 


58.7 


6.6 












0.79 






0.30 






0.91 




60.8 


0.7 




NGC 


2775 


17.03 


(0.02) 


0.96 


1.85 


(0.02) 


0.10 (0.01) 


16.77 


(0.02) 


1.01 


4.60 (0.03) 


37.1 ( 0.4) 


1.4 


3 






16.99 


(0.09) 


0.95 


1.92 


(0.18) 




16.64 


(0.14) 


1.08 


4.40 (0.44) 


29.5 


6.1 












1.07 






0.10 






0.90 




42.4 


1.0 




NGC 


2841 


15.65 


(0.02) 


1.07 


0.86 


(0.01) 


0.47 (0.01) 


16.73 


(0.02) 


1.04 


4.91 (0.04) 


59.8 ( 0.2) 


7.8 


2 






16.19 


(0.08) 


1.06 


1.61 


(0.10) 




17.40 


(0.12) 


1.03 


6.67 (0.55) 


63.9 


6.3 












0.94 






0.30 






1.10 




64.7 


1.2 




NGC 


3593 


15.12 


(0.03) 


1.38 


0.40 


(0.01) 


0.26 


16.98 


(0.03) 


0.87 


2.11 (0.02) 


69.2 ( 0.2) 


10.4 


1 






18.83 


(0.82) 


1.33 


10.16 (22.19) 










3.39 (17.99) 


70.1 


1032.0 












0.86 






0.30 






0.95 




69.5 


0.8 




NGC 


3898 


15.33 


(0.02) 


1.10 


0.90 


(0.01) 


0.22 


17.34 


(0.03) 


0.81 


4.10 (0.06) 


59.7 ( 0.4) 


5.1 


2 






16.50 


(0.15) 


1.09 


2.21 


(0.28) 




18.27 


(0.70) 


0.65 


4.37 (1.75) 


58.0 


7.3 












1.03 






0.25 






0.82 




58.5 


0.9 




NGC 


4378 


17.09 


(0.03) 


0.86 


2.22 


(0.04) 


0.18 (0.01) 


17.62 


(0.04) 


0.88 


4.94 (0.08) 


35.5 ( 0.4) 


1.2 


3 






16.64 


(0.12) 


0.79 


1.60 


(0.18) 




16.72 


(0.17) 


1.06 


3.56 (0.31) 


36.9 


3.7 












0.92 






0.15 






0.78 




37.9 


1.8 




NGC 


4419 


16.58 


(0.03) 


1.36 


1.01 


(0.03) 


0.47 (0.01) 


16.85 


(0.03) 


0.90 


2.40 (0.02) 


73.7 ( 0.2) 


2.2 


3 






17.14 


(0.31) 


1.25 


2.52 


(0.95) 




17.20 


(0.59) 


0.90 


2.44 (0.33) 


70.1 


25.0 












1.19 






0.25 






1.01 




72.1 


1.2 




NGC 


4450 


17.19 


(0.02) 


0.95 


1.44 


(0.02) 


0.25 (0.01) 


17.32 


(0.02) 


0.91 


4.67 (0.04) 


54.6 ( 0.3) 


2.7 


3 






16.77 


(0.10) 


0.83 


1.07 


(0.10) 




16.89 


(0.08) 


1.00 


3.86 (0.24) 


55.2 


7.9 












1.24 






0.05 






0.84 




57.6 


2.7 





Table 2. (continued) 



Name 




(J - K) e 


r 




eb 




(J - K) d 




i 


9 

XI 


n 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


NGC 


4698 


15.91 


(0.03) 


0.76 


0.58 


(0.01) 


0.24 (0.01) 


16.69 (0.02) 


0.96 


3.01 (0.04) 


41.5 ( 0.4) 


7.2 


2 






16.85 


(0.17) 


0.89 


1.12 


(0.17) 




17.69 (0.20) 


0.98 


4.59 (0.72) 


56.6 


15.7 












0.79 






0.08 




0.99 




54.5 


1.8 




NGC 


4845 


16.50 


(0.07) 


1.54 


0.86 


(0.05) 


0.23 (0.02) 


18.05 (0.04) 


0.89 


5.26 (0.10) 


75.7 ( 0.3) 


6.1 


4 






15.85 


(0.21) 


1.59 


0.80 


(0.11) 




18.25 (0.12) 


0.85 


6.03 (0.73) 


75.5 


17.1 












1.43 






0.40 




1.12 




75.9 


1.3 




NGC 


5879 


18.92 


(0.05) 


0.81 


2.45 


(0.05) 


0.18 


16.81 (0.04) 


0.87 


1.40 (0.02) 


75.8 ( 0.3) 


3.3 


3 






16.21 


(1.07) 


0.94 


0.22 


(0.16) 




15.94 (0.09) 


0.85 


1.24 (0.06) 


66.4 


10.0 












1.13 






0.25 




0.74 




65.7 


0.9 




NGC 


6314 


15.79 


(0.03) 


1.16 


1.80 


(0.01) 


0.04 


18.00 (0.06) 


0.83 


8.50 (0.01) 


62.3 ( 0.0) 


2.5 


2 






18.55 


(0.40) 


1.29 


12.15 


(5.06) 








10.62 (18.42) 


70.0 


12.7 












1.04 






0.10 




0.90 




60.9 


1.1 




IC 


724 


16.86 


(0.03) 


1.40 


3.19 


(0.10) 


0.08 


18.46 (0.06) 


0.94 


11.44 (0.44) 


73.1 ( 0.7) 


2.0 


2 






19.54 


(0.33) 


1.43 


23.87 


(9.73) 








12.21 (6.47) 


67.7 


7.7 












1.12 






0.08 




1.14 




70.5 


1.4 





Notes to Table 2. 



Line 1: 
Line 2: 
Line 3: 



2-D parametric decompositions 
1-D parametric decompositions 
if-band 2-D non-parametric decompositions 



BULGE 



DISK 




Table 3. Decomposition results: K-band integrated luminosities and 

bulge-to-disk ratios 



Name 


M b 


M d 


M g 
(4) 


B/D (K) 


B/D (J) 


(1) 


(2) 


(3) 


(5) 


(6) 



NGC 


1024 


-24 


,98 


-24 


.29 


-25 


.44 






1 


.89 






-25 


,17 


-23 


,83 


-25 


.45 






3 


.43 


NGC 


2639 


-25. 


.16 


-25. 


,09 


-25 


.88 


1 


.21 


1. 


.06 






-23. 


,84 


-25. 


,56 


-25 


,76 





.23 


0. 


.20 


NGC 


2775 


-24 


,12 


-25 


.12 


-25 


.48 





.41 





.40 






-24. 


.37 


-24. 


.78 


-25 


.35 





.59 


0. 


.68 


NGC 


2841 


-23. 


.63 


-25. 


,30 


-25 


.51 





.21 


0. 


.21 






-23 


,73 


-24 


.94 


-25 


.25 





.38 


0. 


.33 


NGC 


3593 


-22. 


,16 


-23. 


.21 


-23 


.56 





.24 


0. 


,38 






-21. 


.53 


-23. 


.47 


-23 


.64 





.18 


0. 


.17 


NGC 


3898 


-24 


,04 


-24 


,30 


-24 


.93 





.60 


0. 


.79 






-23 


,77 


-24 


,43 


-24 


.90 





.45 





.54 


NGC 


4378 


-24. 


.46 


-24. 


,41 


-25 


,19 


1 


.06 


1 


.04 






-23 


,91 


-24 


.65 


-25 


,10 





.45 





.50 


NGC 


4419 


-23 


.27 


-23. 


.62 


-24 


.21 





.47 


0. 


.72 






-22. 


,49 


-23. 


,92 


-24 


.18 





.23 


0. 


.27 


NGC 


4450 


-23 


.40 


-24 


.59 


-24 


.91 





.32 


0. 


.34 






-22. 


.90 


-24 


.45 


-24 


.69 





.17 


0. 


.24 


NGC 


4698 


-22. 


.50 


-24. 


.27 


-24 


.46 





.24 


0. 


,20 






-23. 


.22 


-23. 


,94 


-24 


.39 





.62 


0. 


,51 


NGC 


4845 


-23. 


,15 


-24 


,12 


-24 


.49 





.23 





,41 






-23 


,08 


-24 


.03 


-24 


.41 





.31 





,42 


NGC 


5879 


-22. 


,83 


-22. 


,49 


-23 


.42 


1 


.44 


1 


,37 






-20 


,72 


-23 


.25 


-23 


.35 





.07 





.10 


NGC 


6314 


-25. 


,09 


-25. 


,21 


-25 


.91 





.65 


0. 


.90 






-24. 


.89 


-25. 


.39 


-25 


.92 





.55 


0. 


.64 


IC 


724 


-25 


.26 


-25 


.40 


-26 


.09 





.59 





.88 






-25 


.05 


-25 


.58 


-26 


.10 





.63 





.61 



Notes to Table 3. 



Line 1: 2-D parametric decompositions 
Line 2: 2-D non-parametric decompositions 



Bulge parameters 



Disk parameters 
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16 
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Log r e (kpc) Log r d (kpc) 
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Log r (kpc) Log r (kpc) 
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